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Abstract 

Active Galactic Nuclei (AGN) host powerful jets containing high-energy electrons and protons. 
The astrophysical environment where AGNs and their jets are found is characterized by large 
concentrations of both dark matter (DM) and intra-cluster medium (ICM) gas. As the high-energy 
jet particles transverse the DM and the ICM, elastic and inelastic scattering processes generically 
lead to the production of final-state photons. As first envisioned by Bloom and Wells (1998), 
and as more recently pointed out by us and others, the scattering of electrons off of DM could 
lead to a potentially detectable gamma-ray signal, with the parton-level contribution from protons 
offering dimmer perspectives. Recently, Chang et al. argued that taking into account photons 
from hadronization and showering, the actual photon flux is substantially increased. Here, we 
point out that the proton-jets have to be highly collimated, contrary to what predicted by simple 
blob-geometry jet-models sometimes employed in these studies, otherwise they would produce a 
very large flux of photons from inelastic collisions with ICM nucleons, which would outshine by 
many orders of magnitude the signal from DM, for almost any reasonable ICM and DM density 
profiles. 
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Active Galactic Nuclei (AGN) are sources of powerful collimated jets containing high- 
energy particles possibly including electrons and protons, and are believed to exist in en- 
vironments typically hosting large densities of Dark Matter (DM), such as the center of 
massive elliptical galaxies. Ref. [jj pointed out a long time ago that electron-DM scattering 
near AGNs could produce multi-GeV gamma rays, estimated the process from a particle 
physics standpoint and the overall expected photon rate, and concluded that the signal was 
too suppressed to be detectable by current or planned gamma-ray telescopes. 



More recently, Ref. [2| and [3j have elaborated on the original idea proposed by Bloom 
and Wells and have evaluated in detail the scattering of relativistic electrons in AGN 
jets off of DM, arguing that indeed the process can produce a gamma-ray signal potentially 
detectable by the Fermi Large Area Telescope (LAT). From a theoretical and observational 
standpoint, AGN jets can be either of predominantly leptonic or hadronic nature, or they 
can contain both electrons and protons . In the case of hadronic jets, the scattering off 
of DM particles happens at the parton level, and the resulting gamma-ray signal is a few 
orders of magnitude weaker compared to the electron case, after integrating on the parton 
distribution function (PDF). The analyses of both Ref. Q and fl, however, only included 
the parton-level events, and did not include the effects of showering and hadronization. 

In Ref. we argued, based on quantitative estimates, that it appeared unlikely to get 
any significant enhancement from photons produced in the showering of proton remnants. 
Recently, however, Chang et al. reached a different conclusion computing the process in 
detail and including hadronization, and claimed that the resulting gamma ray signal can be 
significantly affected and increased, under certain assumptions on the AGN jet structure. 
The photons from proton-DM collisions are classified as 

1. Final State Radiation (FSR). They dominate at large scattering angles and for photon 
energies up to 5M, typically between a few GeV and ~ 100 GeV, with SM the mass 
difference between the mediator 1 of the scattering process and the DM. Compared 
to the parton-level process, these gamma rays are enhanced by hadronization due to 
large multiplicity and lack of suppression by the QED fine structure constant. 

2. Shower from proton remnants. They are emitted along the proton's incoming direction 



1 In supcrsymmctric models, for example, the DM would be the neutralino and the mediator a squark. 
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and would affect significantly the spectrum at higher photon energies (E 1 > SM). One 
therefore needs a large enough flux of protons in the direction of our line of sight for 
this contribution to be relevant. 

In this short note we point out that if there is a significant flux of protons in the jet 
pointed along our line of sight, the photons produced by inelastic collisions of the jet protons 
with nuclei in the intra-cluster medium 2 (ICM) would by far outshine those produced by 
scattering off of dark matter. We provide here both simple estimates and a detailed worked- 
out example for the case of the nearby AGN Centaurus A (CenA) to substantiate our claim. 

Let us start with a simple estimate, and consider the two processes of proton-DM and 
of proton-proton scattering. The photon flux from the two scattering processes can be cast, 
respectively, as 
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where the label B in the first line stands for "baryons" , while DM in the second line stands 
for dark matter. d^GN is the distance to the AGN ((Lagn — 3-7 Mpc for the case of CenA), 
^b(dm) = / PB(DM)( r ) dr, with Pb{dm)(t) the baryonic (DM) density profile (the factor 
1/mB, dm turns this into the actual number density). da vv( E v) indicates the differential 
cross section for the process pp — > tt + X — > 77 + X, while ^ggMijjgj i s f or the process 



p + DM — y 7 + The last term in the integral, 



, is the energy distribution 

6=0AGN 



of the protons in the jet, in the direction of our line of sight (for CenA, 6agn = 68°). We 
will assume here a simple blob geometry to describe the jet, which implies that the jet is 
not perfectly collimated, but that there are some high energy protons that point away from 
the jet axis. Note that Refs. Q and j^|, while assuming the blob geometry, impose by hand 
a cutoff on the angle of the protons in the black hole frame in order to obtain a collimated 
jet. The authors of Ref. j^], instead, contemplate both a cut-off and a full blob geometry. 
The latter includes "line-of-sight" protons which are traveling towards us and that would be 
responsible for most of the gamma-ray signal from proton-DM scattering at energies higher 
than 5M, as we mentioned above. 



2 While the ICM consists primarily of ionized hydrogen and helium, we shall assume here for simplicity 
that it be composed of ionized hydrogen (i.e. protons) only. 
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Let us first compare the relevant differential photon fluxes for the two processes (i.e. from 
scattering of high-energy cosmic-ray protons from the jet off of ICM protons versus off of 
dark matter). The overall factors appearing in Eq. ([T]) for proton-DM versus proton-ICM 
proton differ by three elements: 

(i) for proton-DM scattering, the total cross section is at most 10 -5 mb 5], while for the 
inelastic process pp — > ir + X, the inclusive inelastic cross section is about 30 mb, i.e. six 
orders of magnitude larger. 

(ii) since the 7-ray flux depends on the number density of targets, and not on the physical 
density, the fluxes are inversely proportional to the mass of the target particles, which for 
DM masses of order 100 GeV, as contemplated in these studies, gives an extra suppression 
of two orders of magnitude for proton-DM compared to proton-proton. 

(hi) the last, critical element is the comparison of the DM and baryon density profiles. 
While there exist large observational uncertainties, numerical simulations indicate that the 
physical densities of baryons and DM are, over the relevant radii, comparable 3 . Assuming 
that is the case, the values of 8b and 5dm are of the same order. We will quantitatively test 
caveats to this assumption (which might depend on a steep density profile for r — > 0) below. 

Thus, from this preliminary estimate, we expect the photon flux from proton-proton 
scattering to be several orders of magnitude (between 6 and 8, but possibly more) larger 
than the one from proton-DM. If indeed one allowed a jet geometry with a fraction of off-axis 
protons, as in the naive blob model without a cutoff, one would conclude that the photon 
flux from proton-ICM proton inelastic scattering would completely outshine the process 
involving dark matter and violate gamma-ray observations by many orders of magnitude. 

We intend to test here the robustness of the conclusion we reached in the qualitative 
discussion above, and assess possible caveats to that conclusion by contemplating the pos- 
sibility that a very steep DM density profile compensate for the many orders of magnitude 
larger 7-ray flux expected from scattering off of gas nuclei. We parametrize the density pro- 
files for the baryonic and dark matter target densities for small radii r — > by the following 



3 See for example figure 4 of (6|. To test this in the innermost regions (where the jet structure becomes, 
however, increasingly uncertain) one would have to extrapolate to much smaller radii than what showed 



e.g. in Ref. [6| , but the trend indicates that the density of cold gas is generally higher than that of DM 
as we approach the center 
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simplified functional forms: 



r 



Pdm — Po,dm — , (3) 
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While we assume the same scaling radius, 5b, dm will simply scale linearly with an overall 
normalization change induced by a change in r c for either one of the two species, so taking this 
effect into account is trivial. The key aspect is what happens when a and (3 are drastically 
different. 

We normalize the overall factors po r " ,/? to be of the same order of magnitude as the 
mass enclosed within 5 kpc from the AGN: in the case of CenA the total gravitating mass 
enclosed in 5 kpc is equal to 7.5 x 10 10 M & , as obtained from X-ray measurements under the 
hypothesis of hydrostatic equilibrium Notice that this is, for our purposes, a conservative 
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assumption, as the DM density is generally thought to be subdominant (i: 
compared to the baryonic mass in the vicinity of the AGN center (see e.g. 
the factor 5b in Eq. ([!]), we set r min = 10 Rs, with Rs = 5 x 10~ 6 pc the Schwarzschild 
radius for CenA, and r max = 15 kpc. The integral has a very mild dependence on r m j n when 
the latter is chosen between a few and 1000 Rs, while it is very insensitive on the upper 
limit r max (which depends on how far the jet extends from the AGN). 

In the left panel of Fig. [1] we vary the steepness of the density profiles a and (3 for the 
dark matter (a) and for the baryonic (/3) densities, and show curves of constant 5b,dm- The 
labels indicate curves at the corresponding values of log w (5r)M) (white labels) and log 10 (^s) 
(red labels). The contours saturate at close to constant values towards the lower left corner. 
As we anticipated in the qualitative discussion above, for comparable steepness a, (3, we in 
fact find 5dm ~ 5b- 

We now proceed to the calculation of the 7-ray flux from jet protons scattering off of DM 
versus ICM protons. For the DM calculation, we closely follow the assumptions in Ref. j^J, 
and we employ here the following fiducial values for the dark matter particle properties: a 
mass (mfljtf) of 300 GeV, and a total cross section (<7 p dm) of 10~ 5 mb. 

To calculate the photon flux from proton-proton interactions, the relevant differential 
cross section can be written as [4j 

d ^^^2a pp ,^{E p )5{E 1 - X E P l (5) 
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FIG. 1: In these figures a and (5 are the slopes of the DM and baryon profiles, respec- 
tively. In the plot on the left the vertical lines are contours of log 10 (<5DM)) while the hor- 
izontal lines are contours of log 10 (5_B)- On the right we show contour lines (thick solid) of 
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10 GeV. The two vertical dashed lines show that 
can be below the observational limits for lower val- 
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ues of the slope a. The red line corresponds to a flux of 5 x 10 -12 erg s _1 cm~ 2 , value measured 
by Fermi at 10 GeV, and the shaded region above it is then excluded. 

where 0" PPi7ro ~ 30 mbarn is the inclusive inelastic cross section for the process pp — > tcq for 
energies E p > 10 GeV. The delta function is an approximation for gamma rays formed in 
the reaction pp — )> 7r — )> 77, with \ ~ 0.05 the ratio of the mean energy of the produced 
photon to the energy of the incident proton. 

Assuming a simple blob geometry to describe the cosmic-ray proton jet, and again em- 
ploying the same assumptions as in Ref. j^], we write the proton energy distribution in the 
black hole frame as 
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where r# = 3 is the bulk Lorentz factor (of the blob with respect to the black hole), which 
corresponds to a velocity 0b = 0.94, L p = 10 45 erg/s is the luminosity of the jet, E min = 10 



GeV and E n 



10 7 GeV are the minimum and maximum energies of the protons in the 



distribution. E Q is some reference energy that does not need to be specified, as it cancels 
out in the definition above. As a result, we can cast Eq. ([T]) in the following form: 

We note that the quantity typically shown is the spectral energy distribution, i.e. vS v = 
E 2 ^-: here, vS v is close to a constant, as the powers of E 1 cancel out. 

The right panel of Fig. [T] illustrates curves of constant 7-ray flux at an energy E~ = 10 
GeV. We shade the region ruled out by Fermi observations of CenA. We note that the entire 
plot features fluxes in excess of 10 -12 erg s -1 cm -2 due to the proton-proton contribution. 
As can be appreciated from the dashed vertical lines, for the DM contribution to overcome 
the smallest possible proton-proton background one would need DM density profiles much 
steeper than a = 2, a value that is extraordinarily hard to achieve in realistic situations, 
especially in the (most relevant) innermost regions: there, tidal disruption due to stars 
would very likely make an extremely steep density profile highly unstable. Also, the very 
large densities achieved at small radii for steep profiles with a > 2 would likely be cutoff 
by DM pair annihilation [9( at a radius potentially much larger than the value of R m in than 
the one we assumed to calculate 5dm here. 

Besides the steepness of the density profiles, an additional assumption one might consider 
changing is the picture for the jet structure and geometry. Remember that with the blob 
geometry, although most of the protons move within a cone along the direction of the jet- 
axis, there is still a small fraction pointing at large angles with respect to the axis. This 
small fraction leads to the catastrophic result just described, so one might wish to suppress 
it. Indeed, observations indicate that AGN jets appear to be very collimated. One way to 
suppress off-jet protons, is to introduce "by hand" a cut-off on the angle in order to restrict 
the particles to move in a narrow cone centered on the axis, as was done e.g. in [2|. In 
Figure 7 of Ref. |5| the authors show results with such a cutoff (orange lines) and without 
it (black line). In light of our considerations, we argue that the latter should be discarded. 

As a further remark, we wish to point out that if the jet is well collimated, photons 
from proton remnants should be negligible. Let us consider the production of energetic 
7r°'s at large angles, e.g. 9 « 68°, in the case of a collimated proton jet scattering off 
ICM protons. This process requires one active quark to be bent at a large angle and it 
involves a 4-momentum transfer in an elementary collision of order t ~ — E 2 , where E 



denotes the scattered quark energy. The amplitude for such a process would be suppressed 
as (A/E) 2 , A ~ 1 GeV being a characteristic hadronic scale; the differential cross section, 
correspondingly, would be suppressed as da/dVL{6 = 68°) ~ da /dVt{6 = 0)(A/E) 4 , implying 
that for E y ~ 100 GeV we would encounter a 10~ 8 suppression relative to n° production in 
the direction of the AGN jet. 

Note that for a predominantly leptonic jet, i.e. the case that Ref. and Q argued being 
the most promising one, the relevant associated background of jet electrons scattering off of 
ICM protons producing bremsstrahlung photons would be by far suppressed compared to 



the proton-proton background we studied here 



.a 



101 ] . The present discussion, therefore, does 



not invalidate the results of Ref. 1211 and |3[ for a predominantly leptonic jet. 

In conclusion, we studied the impact of inelastic proton-proton collisions in AGN jets that 
could produce a sizable gamma-ray signal from proton-DM scattering. While Ref. [5[ argued 
that photons from proton remnant showering could contribute to enhancing the proton-DM 
signal, particularly at higher photon energies, we showed here that those same jet protons 
would vastly outshine the proton-DM gamma-ray flux as they inelastically scatter off of 
nucleons in the intra-cluster medium. We argued that the proton-proton background is at 
least 6-8 orders of magnitude larger than the DM-proton signal for comparable baryonic and 
dark matter density profiles. 

Finally, while the process of dark matter scattering off AGN jet particles is potentially 



observable if the AGN jet is predominantly leptonic (as pointed out in Ref. |2j and [3]), we 
remain somewhat pessimistic as to the possibility of detecting a signal for the DM-proton 
case. 
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